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ABSTRACT 



Discrete symmetries in grand canonical ensembles and in ensembles canon- 
ical with respect to triality are investigated. We speculate about the general 
phase structure of finite temperature gauge theories with discrete Z(N) sym- 
metry. Low and high temperature phases turn out to be different in both 
ensembles even for infinite systems. It is argued that gauge theories with 
matter fields in the fundamental representation should be treated in ensem- 
bles canonical with respect to triality if one wants to avoid unphysical pre- 
dictions. Further, we discuss as a physical consequence of such a treatment 
the impossibility of the existence of metastable phases in the quark-gluon 
plasma. 
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The recent interest in studying the Z(3) structure of hot gauge theories 
is caused by the strong will to obtain more insight into the phase 
structure of strongly interacting matter for physical models and for specu- 
lations. In the present article we follow the obvious tendency of the above 
papers to avoid a polemic with each other, to supplement the conventional 
approach |10| and to develop our own point of view |], |2[. The guideline 
of this exploration is to draw the arguments supporting the conclusion of 
[P], |3j which come from the impressing recent activity on the role of local and 



global discrete symmetries in gauge theories fll], [L^, [L3], [14]]. Usually, these 
symmetries appear as a commuting centre of an underlying gauge group. An 
important conclusion of the above papers is an appearance of a new Higgs 
phase where the corresponding centre charges are not screened while a gauge 
symmetry is broken up to a discrete subgroup. Moreover, unlike in the con- 
fined phase these charges could, in principle, be detected at long distances 
via the Aharonov-Bohm effect. In order to make our arguments more clear 
we have first to review some aspects of the present status of MC analysis of 
finite temperature QCD. 

1 Z(3) symmetry in the grand canonical en- 
semble description (GCE). Euclidean for- 
mulation 

In lattice QCD the gluon fields appear in the form of SU (3) matrices U x ^ 
which are defined on links (x, fi) of a four- dimensional euclidean lattice. In 
finite temperature lattice QCD the correlation function < L{fi) ■ ■ ■ L*{f N ) > 
of several Polyakov loops 

i N t 

L(r) = -Trl[U (r,t) (1) 

6 t=i 



can be connected |10]] with the free energy F of iV infinitely heavy quarks q 



or antiquarks q at the corresponding positions r*i, • • • , r N at temperature T 
relative to the free energy of the vacuum 

F(q(n), - ■ ■ ,q(r N )) = -T\n < L(n) ■ ■ ■ L*(r N ) > . (2) 
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The thermodynamical average < • ■ ■ > is computed using the partition func- 
tion 



Z = / V[U] e 



-S[U] 



(3) 



as a "sum" over all gauge field configurations U. In pure gluonic QCD for 
most problems the appropriate choice for S is the Wilson action 



where the plaquettes U x ^ v are built of four links in the /iz/-plane of a four 
dimensional euclidean lattice. 

The action Sa[U] is Z{3) symmetric. This means that multiplication of 
all links in direction ji = in the three dimensional x, y, z-torus with fixed 
t, e.g. t = by a Z(3) element leaves the action invariant. Such a global 
Z(3) transformation can be regarded as an aperiodic transformation of the 
centre of the gauge group. A single Polyakov loop L(r) transforms under 
Z(3) nontrivially and therefore the distribution of L(r) values can be used 
as an indicator for spontaneous breaking of Z{3) symmetry. 

Let us now discuss the results of lattice calculations and their common 
interpretation. Monte-Carlo calculations in pure gauge lattice QCD show 
usually a characteristic behaviour of the spatial average L of L(r) during 
a Monte-Carlo simulation. In the confinement regime (low j3) L scatters 
symmetrically around zero in the complex plane. In the deconfined phase 
there appear three Z{3) symmetric maxima of the Polyakov loop distribution 
in 0°- and ±120°-directions. The appearance of the three peak structure of 
the L-distribution in the deconfinement regime close to the phase transition 
is considered as a demonstration that spontaneous symmetry breaking on a 
finite lattice can never happen exactly. As the number of tunnelling events 
between the maxima decreases with increasing (3 one commonly expects L 
to be frozen in the thermodynamical limit in one of the Z(3) directions and 
thus to get spontaneous Z(3) symmetry breaking. Therefore, one may obtain 
arg < L >= 1, ±^ in the thermodynamical limit. With dynamical fermions 
the action S contains a fermionic contribution Sf 




(4) 




x,^ x+p,,v u x+U,fl U x,v ' 
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S — Sq + Sf, 



x,x' 



which breaks Z(3) symmetry explicitly. In the Kogut-Susskind formulation 
[ |15|j the fermionic matrix M reads 



Mt,i' — 9 (r , ! r,/it4,/i^a:',a;+/u — ^ x' ,pUl' ^x' + Tn5 x>x r, 



p _ f-nzi+a^H hx M -i 



(6) 



with one-component Grassmann variables ty x and ty x . The integration over 
the Grassmann variables in the path integral 



Z = I V[U,V,V] e~ s (7) 
can be performed analytically and leads to the fermionic determinant detM 

Z = J V[U] e- SG[u] detM. (8) 

It is easy to see that the fermionic action violates Z(3) symmetry [|l^ . 
In the usual MC-iterations one finds that in the low temperature phase the 
expectation value < L > of a Polyakov loop L is a small positive number 
and in the high temperature phase the three peak structure of the distri- 
bution becomes asymmetric. The maximum at 0° is favoured and the two 
lower maxima at angles of around ±120° are symmetric. This fact is usually 
interpreted as a consequence of the violation of Z(3) symmetry in QCD. 

Thus, the lore is: 

In the low temperature phase Z(3) symmetry signals the confinement of 
colour charges and in the high temperature phase spontaneous breakdown 
of Z{3) symmetry means colour charges screened by gluon fields and the 
Polyakov line is an appropriate order parameter. In QCD with dynamical 
quarks, Z{3) symmetry is explicitly broken and Polyakov loop ceases to be 
the appropriate order parameter. It should be stressed that these results 
reflect the grand canonical ensemble treatment of full QCD with respect to 
triality. 
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These results are grounds for nonperturbative MC analysis of QCD pro- 
viding a theoretical basis of quark-gluon plasma phenomenology. However, 
this treatment is not complete and absolutely satisfactory. There exist sev- 
eral reasons to believe in such a conclusion. This approach does not solve 
the long standing problem of an order parameter for full QCD. Moreover, 
it fails in tempting to treat Z(3) metastable phases as physical states in 
Minkowski space due to running into thermodynamical inconsistency. We 
shall argue that this picture is conceptually not perfect as well, since it does 
not substantiate unified mechanism of quark confinement at zero and finite 
temperature. Indeed, confinement at zero temperature implies only zero tri- 
ality states and unbroken Z{3) symmetry which triality is connected with. 
Meanwhile, the corresponding symmetry is explicitly broken in full QCD at 
any low temperature. At least, it does not answer the important question 
what occurs if the triality charge is not screened and develops - according 
to the appropriate Gauss law - a long-range observable which could be, in 
principle, experimentally measurable. 

We are going to demonstrate below that there are no such problems in 
the canonical ensemble description of finite temperature QCD. 

2 Z(3) symmetry in the canonical ensemble 
description (CE). Euclidean formulation 

We develop a view Jl|, [|, [|] different from that mentioned above which seems 
more consistent at least in many general aspects. Let us begin mentioning the 
Euclidean formulation of the theory at zero temperature. The Lagrangian of 
QCD is symmetric under transformations of the gauge colour group SU(3)i 0C . 
There is a global subgroup of the local symmetry which can be decomposed 
as SU(3) g i = SU{3)/Z(3) g i ® Z{3) g i. Gauge fields are identically invariant 
under Z(3) g i transformations but fermionic fields transform as 

-> e i2 r k ^(x) = e^ A8fc ^(x) (9) 

leaving the QCD Lagrangian invariant. This symmetry is related to triality 
t q which is defined as 

t q = (N q - N q )mod 3 (10) 
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or as expectation value of the operator 



t q =< ^(x)V3\ 8 ty(x) > . (11) 



It had been proposed [|17| that in QCD the gluon colour charges and any zero 
triality charges are screened by a dynamical Higgs mechanism (understood as 
a screening mechanism) leaving Z(3) symmetry unbroken and a single quark 
(having nonzero triality) unscreened. Then long-range chromoelectric forces 
between quarks will persist according to the Gauss theorem for a discrete 
centre. 

In the scheme described in the previous section (GCE) a single quark may 
have a finite free energy and the states with open triality can contribute to the 
partition function. As a result the domains with different Z(3) phases give 
different contributions to the partition function. The zero phase contributes 
most. We are sticking here to the opinion that this problem is an artifact of 
the GCE description which leads to the explicit violation of the Z(3) ape- 
riodic symmetry. We shall show below that these aperiodic transformations 
at finite temperature can be regarded as acting on links in time direction or 
equivalently on fermionic fields. These finite temperature transformations of 
fermionic fields must be treated in analogy to the zero temperature transfor- 
mations (^]). They are related to triality in the same manner as (Q). Since 
our premise is the statement that QCD with dynamical fermions is able to 
describe only systems of zero triality as it is at zero temperature and all 
quark phases should be indistinguishable both in the partition function and 
in thermodynamical functions we would like to construct a theory which is 
invariant under the corresponding Z(3) transformations. 

Let us first describe how one can construct such a theory for pure gluo- 
dynamics. In the absence of external perturbation (like "magnetic" field for 
spin systems), the average value of the Polyakov loop will be equal to zero 
at any temperature even in the thermodynamical limit: as J T>[U] e~ Sa ^ is 
Z(3) invariant and the Polyakov loop ([I]) is not, one has to get 

<L>=0 (12) 

in the confined and also in the deconfined phase. The important point is there 
is no contradiction between this statement and the spontaneous breaking of 
Z{3) symmetry in the deconfined phase. To demonstrate that one can use 
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< L(0)L*(f) > correlations or even more suitable the {Z(3) symmetric) L(f) 
distributions expressed by the moments |L(r)| 2 , |L(f^)| 3 , • • •, what is equiva- 
lent to the determination of Polyakov loops in Z(3) invariant representations 
(octet, decuplet, antidecuplet, etc) fl8fl . 

An obstacle for accepting < L >= in the deconfmed phase comes from 
lattice calculations which give < L >^ at high (3. This result has to 
be compared with the usual method to demonstrate spontaneous symmetry 
breaking by an external field. Including —A J2r Re{zL(f)}, z G {1, e ±2nl ^ 3 } in 
the QCD action, performing the thermodynamic limit and then the limit to 
vanishing A would result in < zL >^ and positive in the deconfmed phase. 
As we cannot influence a physical system by a field which we have switched 
off we have to give an interpretation of this result: The external field fixes 
the coordinate system in which the "direction" of spontaneous symmetry 
breaking is measured. This direction is unchanged even after switching off 
the external field. To get < L >= we have to use the "unbiased" path 
integral which averages over all directions of such a coordinate system and 
gives < L >= also in the deconfmed phase. The observation that for 
a Z(3) invariant Lagrangian the observables of nonvanishing triality have 
expectation value zero is rather trivial for the pure gluonic case (see Jl], || for 
more details). Nevertheless, the usage of the "unbiased" path integral (i.e., 
with perturbation and averaging as described above) is more preferable since 
we do not have violation of triality in such a treatment despite deconfinement 
takes place. We will now turn to the more interesting case of full QCD. 

Periodic boundary conditions in imaginary time direction imply that we 
are dealing with a theory defined in a space with topology R 3 <g> S 1 . To con- 
struct a theory with conserved triality and indistinguishable phases of quark 
fields we will follow a recently developed ideology Q, |13| (see also [0, [L4| 
and references therein). Fields of gauge theories on space-time manifolds that 
are not simply connected need not to be strictly periodic on noncontractible 
closed loops of this manifold. Instead, these fields need only be periodic up 
to the action of an element of Z{3). Just these boundary conditions distin- 
guish local Z(3) symmetry from global Z(3) symmetry and allow the local 
symmetry to manifest itself through a nontrivial Aharonov-Bohm effect. 

Let us consider the Taylor expansion of exp(— Sp) in the path integral. 
The first nontrivial term, coming from time-like links, which will survive the 
integration over fermionic fields is proportional to fit ^t^tUo(t). This term 
describes propagation of a single quark along the temperature direction. The 
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triality transformations of the fermionic fields in ]J t ^f t ^/ t Uo(t) should give a 
nontrivial phase factor. It is obvious that one can achieve this by considering 
in addition to (|9[) the following substitutions for the fermionic fields 

-> -e i2 -? k V{x,0). (13) 

If we now perform this substitution in the operator Y\ t ^t^tUo(t) and then 
consider the action of the transformation @ on this operator we will find 
that the operator gets a nontrivial phase factor as it should be for a state 
with a single quark. Thereby, this is the obvious extension of triality trans- 
formations (0) to the finite temperature theory. The phase transformation 
(H) determines the triality of quark states at zero temperature. Hence, this 
transformation together with (|T3"| ) determines the triality of quark states at 
finite temperature. 

In the zero temperature theory the transformation (|9|) acts trivially on 
all physical states. We would like to construct a formulation of QCD where 
the corresponding finite temperature transformations act trivially on physical 
states as well. To get such a theory where only triality zero states contribute 
we have to sum over the following boundary conditions in the time direction 

V(£,/3) = -e^ fc ^(f,0). (14) 

This is equivalent to a summation over all phases in (|9]) combined with (|i~3|) 
and to a projection onto zero triality states. 

We explain now that such a projection is nothing but the requirement of 
local gauge invariance in the meaning of refs. and [O. The boundary 



conditions flT4]) allow to identify field configurations which are related by the 
transformation (liHT). Such a relation can be considered as the definition of a 



local discrete symmetry [I2| and is known as the so-called orbifold construc- 
tion. It was shown [12] that one must project out noninvariant states from 
the theory with local gauge symmetry of such a type. Thus, because of the 
identification (|i~4D all physical observables should be invariant under trans- 
formation of the local Z(3) group. All noninvariant states, like single quark 
states, must be projected out from the theory to obtain a true Hilbert space. 
This defines an ensemble canonical with respect to triality. On the other 
hand, it is a natural extension of the orbifold construction for nonabelian 
gauge theories at finite temperature. 
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There is no difficulty to demonstrate that the summation over all bound- 



ary conditions (14) is identical with the projection onto the zero triality 
sector of the fermionic determinant (a similar proof that only quark loops of 
zero triality contribute to the partition function in this case has also been 
given in Q). We can define the triality operator T (with eigenvalues 0, ±1) 
which determines the triality of closed loops, and with his help a projection 
operator Pq- on triality T 

P T =- £ e fc2m(t " r)/3 . (15) 

^ fc=0,±l 

The implementation of this projection operator in the lattice formulation is 
very simple. One has to multiply all links in time direction for an arbitrary 
time step with the phase e k2m ^ 3 and to determine with this new link variables 
the determinant. The result has to be multiplied with e k2mT / 3 _ Finally, 
the sum over phases k — 0, ±1 has to be executed. Then the triality zero 
contribution of the fermionic determinant is 

det M = P detM = - ^ e k27rif/3 detM (16) 

^ fc=0,±l 

and the full fermionic determinant is the sum over all triality sectors 

detM= det r M, (17) 

r=o,±i 

with 

det r M = P r detM. (18) 

It can be easily seen that this procedure is exactly equivalent to the sum- 
mation over all boundary conditions (H). From eq. ([D]) it is obvious that 



the observable L having triality T = 1 tests only the T = — 1 sector of the 
fermionic determinant which is not the vacuum sector. Using the correct 
T = vacuum sector for the fermionic determinant results in a Z(3) sym- 
metric L-distribution. Let us discuss this in more detail. As the measure 
T>[U] and the gluonic Lagrangian Sg[U] are Z(3) symmetric, the product 
det M[U}0 T [U] in 



<o r > z 



l - J V[U}e- SG[u] detM[U]Or[U] (19) 
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has to be also Z(3) symmetric to get < Or >^ 0. In other words Z(3) 
violating contributions in detM[[/]Or[^] are automatically eliminated by 
the path integral. Only the triality — T component det_rM of the fermionic 
determinant survives. det_r^f has T quark loops less than antiquark loops 
winding around the lattice in time direction. 

Applying this statement to the example of the single Polyakov loop L 
representing an infinitely heavy quark Q, the fermionic determinant supplies 
a triality —1 state of light quarks which will mostly consist of a light antiquark 
q in order to colour neutralize the heavy quark Q. The Polyakov expectation 
value < L > is a thermodynamical mixture 

< L >= ... + e ~ Fm/T + e - F ^ qq),T + ... (20) 

of a heavy-light meson Qq, a heavy-light baryon Qqq, ... . These heavy-light 
states should be bound in the hadronic phase and "ionized" in the quark- 
gluon plasma phase. The charge density of such light quarks around static 
sources has been measured in ref. . 



For the system of a heavy quark Q at position f\ and a heavy antiquark 
Q at position r 2 

L{n)L*{f 2 ) (21) 

a totaly different component of the fermionic determinant, the triality zero 
component, contributes. It is evident that this component is Z(3) symmetric, 
i.e. the three Z{3) transformed sectors which differ in L by factors e ±_ 5~ 
contribute with the same Boltzmann factor. This example shows that for 
observables of triality zero full QCD is Z{3) symmetric. The fermionic action 
does not destroy Z(3) symmetry if we choose to work in the present scheme. 



3 Z(3) symmetry in the Hamiltonian formu- 
lation 

It may be instructive to see how the Z{3) symmetry enters the Hamiltonian 
formulation of QCD. In temporal gauge we have to integrate over all time- 
independent gauge transformations of gauge and of matter fields. This defines 
a projection operator 

P = f D[a £ ] exp[-z f dxa a (x)((DE) a - p a ) + % f dsa a E a ], (22) 
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where p a = \j/^(A a /2)\I/£. P projects onto "physical" states which have to 
fulfil a local Gauss law: (DE) a = p a and a global one: f v divE a dV = Q a = 
0. The global Gauss law warrants that physical states are global colour 
singlets. We would like to concentrate here on the centre of the global gauge 
transformations. In this case the operator (|22f) takes the form 

P = 5((DE) a - p a )P tr (23) 

which includes a Kronecker-delta for the quark triality charge 

, k=N-l Oaroh 

P* = j; E ex P [^(E^A^)] (24) 

iV fc=0 iV x 

(where A a = \^3\ 8 for SU(3)) and, thus, defines canonical ensemble with 
respect to triality. Following the usual procedure one can return to the 



Lagrangian formulation with a summation over all boundary conditions ([LJ) 
for fermionic fields. 

Instead of the canonical ensemble with B = T = which should be well 
suited, a grand canonical description is allowed if it predicts the same be- 
haviour in the thermodynamical limit. Common lattice Monte-Carlo calcula- 
tions in full QCD use the full fermionic determinant detM and thus simulate 
a grand canonical ensemble with chemical potential p = containing all three 
triality sectors. We want to point out that there are neither requirements nor 
conditions in the theory of strong interactions with long-range gauge fields 
guaranteeing that these two descriptions should coincide pi, 22\. One of the 



main conclusion of ref. ^TJ is that one must verify the agreement between 



ensembles in each particular case (i.e., for each conserved quantum number) 
separately. It has been proved |2l| for free theories that canonical and grand 
canonical ensembles are identical in the thermodynamical limit with respect 
to the baryonic number and the strangeness. A similar proof for the global 



colour charge has been done in ref. |E3l . We have shown here that for the full 



theory grand canonical and canonical ensembles differ with respect to quark 
triality: the Z(3) symmetry has different realizations in these two ensembles. 
In ref.[|J it was demonstrated that some of thermodynamical functions may 
have different behaviour as well. 

The most often objection against the above picture is that the projection 
onto the zero triality sector defines a new theory which can somehow dif- 
fer from the genuine QCD. We believe that the present consideration avoids 
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this misunderstanding, which has nothing to do with the discussed situation. 
We have shown that this projection is equivalent to using the canonical en- 
semble with respect to triality in the Hamiltonian formulation. Further, the 
Lagrangian formulation of the theory showed that we did not modify the 
QCD action. The canonical ensemble description is realized in this formula- 
tion by means of a summation over the corresponding boundary conditions. 

4 Phase structure of full QCD 

The present consideration allows us to reexamine a possible phase structure 
of full QCD and to speculate about the high temperature region within the 
canonical ensemble description. One may claim by now that a quantum field 
system with a nonabelian local gauge symmetry including discrete (local 
and/or global) Z(N) symmetry can be in one of the following states depend- 
ing on the concrete representations of gauge and matter fields in the initial 
Lagrangian: 

1. Coulomb phase, gauge fields are massless. 

2. Higgs phase of the first type, triality screening phase. Gauge symmetry 
is completely broken, gauge fields acquire masses. 

3. Higgs phase of the second type, all colour charges excepting nonzero 
triality are screened. Gauge symmetry is broken up to its discrete centre, no 
massless gauge fields. Triality charges are screened classically but may be 
detected at long distances by specific quantum mechanical process. 

4. Confinement phase. Triality charges are unscreened in any sense but 
cannot be detected at long distances. Gauge symmetry is broken either up 
to its Z(N) subgroup (vortex mechanism of confinement) or up to its abelian 
subgroup (monopole mechanism). 

We would like to conjecture the following scenario to be realized in QCD 
with dynamical quarks in the fundamental representation: 

1. QCD at zero temperature always is in the confining phase. It seems 
this statement is generally accepted by now although there is no strict proof 
of this fact. 

2. Full QCD at finite temperature, as described in the grand canonical 
ensemble, always is in the Higgs phase of the first type where also triality 
charges are screened. The aperiodic Z(3) symmetry is explicitly broken at 
any temperature and domains with different phases contribute in different 
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ways to the partition function. In the thermodynamical limit only one do- 
main with fixed phase (k = in the commonly used description) can be 
found. 

3. Full QCD at finite temperature, as described in the canonical ensemble, 
is in the confining phase at low temperature. When the temperature increases 
a phase transition to the Higgs phase of the second type takes place. To 
decide what phase is realized in both ensembles one can apply an order 
parameter proposed in ||13|| . This order parameter can probe a phase of a 
matter field which carries Z(N) charge. In the GCE only one domain can be 
found in the thermodynamical limit. It implies that the corresponding order 
parameter will be always equal to 1 (k = phase). In the CE all domains 
with different phases are degenerate (see discussion just below) and we should 
expect that this order parameter can distinguish between different phases in 
the thermodynamical limit. This means that triality charges of quarks can 
be detected at long distances via the Aharonov-Bohm effect [II], [T3|j . 

One of the interesting consequences of the canonical description concerns 
the possible states of quark-gluon matter at high temperature. Recently, 



in refs. |24] , |25| it has been concluded on the one hand that explicit Z(3) 
symmetry breaking leads to the existence of metastable states at arbitrary 
high temperatures. On the other hand the authors of refs. [^, [27j be- 



lieve that Z(3) phases cannot be prepared as real macroscopic systems [ p8| 
These models are based on the Z(3) asymmetric L-distribution in the quark- 
grand canonical ensemble. But as we argued above one has to investigate 
the zero-triality sector of the fermionic determinant in order to study the 
QCD vacuum. In the corresponding ensemble one obtains Z{3) symmetric 
L-distributions and hence no metastable phases can be found which persist up 



to infinitely high temperatures. The effective models of refs. [24], |25|, |26| \27j 
are effective models of heavy- light hadrons in the QCD vacuum and not effec- 
tive models of "the" QCD vacuum. Exact two- loop calculations in CE, which 
demonstrate that all Z(N) phases are left degenerate in CE, have recently 
been presented in ||29||. It is just this degeneracy which gives the possibility 



to observe different Z(N) domains via the Aharonov-Bohm effect. 
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